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ABSTRACT 

We test whether the spiral patterns apparent in many large disk galaxies should be thought of as 
dynamical features that are stationary in a co-rotating frame for > tdyn^ as implied by the density 
wave approach for explaining spiral arms. If such spiral arms have enhanced star formation (SF), 
observational tracers for different stages of the SF sequence should show a spatial ordering, from up- 
stream to downstream in the corotating frame: dense H I, CO, tracing molecular hydrogen gas, 24 /im 
emission tracing enshrouded SF and UV emission tracing unobscured young stars. We argue that such 
a spatial ordering should be reflected in the angular cross-correlation (CC, in polar coordinates) using 
all azimuthal positions among pairs of these tracers; the peak of the CC should be offset from zero, in 
different directions inside and outside the corotation radius. Recent spiral SF simulations by Dobbs 
& Pringle, show explicitly that for the case of a stationary spiral arm potential such angular offsets 
between gas and young stars of differing ages should be observable as cross-correlation offsets. We 
calculate the angular cross-correlations for different observational SF sequence tracers in 12 nearby 
spiral galaxies, drawing on a data set with high quality maps of the neutral gas (H I, THINGS), 
molecular gas (CO, HERACLES) along with 24 /im emission {Spitzer, SINGS); we include FUV images 
{GALEX) and 3.6 /im emission (Spitzer, IRAC) for some galaxies, tracing aging stars and longer 
timescales. In none of the resulting tracer cross-correlations for this sample do we find systematic 
angular offsets, which would be expected for a stationary dynamical spiral pattern of well-defined 
pattern speed. This result indicates that spiral density waves in their simplest form are not an 
important aspect of explaining spirals in large disk galaxies. 
Subject headings: galaxies: general — 



1. INTRODUCTION 

The nature of a coherent spiral structure that extends 
over a large portion of the galaxy is still not fully un- 
derstood. While there are many variations in theories to 
explain such a structure, typically self-excited models for 
spiral structures can be divided into two groups based on 
the longevity of the spiral pattern. A long-lived quasi- 
stationary spiral structure theory has been developed 
largely using analytical studies. The spiral features are 
attributed to quasi-steady global modes of the disk (Lin 
& Shu 1964, 1966; Lynden-Beh & Kalnajs 1972; Bertin 
et al. 1989a, 1989b, etc.). The other approach, which has 
largely been based on simulations, purports that spirals 
are short-lived, recurrent, transient patterns that orig- 
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inate from recurrent gravitational instabilities (Goldre- 
ich & Lynden-Bell 1965; Julian & Toomre 1966; Toomre 
1981; Sellwood & Carlberg 1984; Sellwood 2010). It re- 
mains a challenge to decipher which of these descriptions 
best describes spiral structure. The interested reader can 
refer to excellent reviews in Binney & Tremaine (2008) 
and Sellwood (2010) for more information. In these re- 
views it is explained that it has proven very difficult to 
decipher which of these models best describes observed 
galaxies and it may be that a combination of these mod- 
els is required. It is clear that careful comparisons be- 
tween observations and simulations are needed in order 
to disentangle these theories. 

In this work, we examine simulations and observations 
for evidence of a long-lived, quasi-stationary spiral struc- 
ture based on the predictions of such a model. As origi- 
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nally explained by Roberts (1969, hereafter R69), if the 
spiral pattern is a quasi-stationary, dynamical pattern in 
a corotating frame and SF is occurring at an enhanced 
rate in the arms, then there should be a temporal se- 
quence of events, as material streams in and out of this 
spiral pattern. 

As gas passes through the minimum of the spiral poten- 
tial, it gets compressed and molecular clouds form. This 
is followed by SF, leading initially to dust-enshrouded 
young stars and then to unobscured young stars and 
clusters. If this SF sequence of events were to occur 
in a steady state, then it would translate into a set of 
spatial offsets for tracers of different stages of this SF 
sequence, with the size of the offsets reflecting both the 
time difference between these stages and the relative ve- 
locity between the material and the spiral pattern. In 
this way, observed spiral arms are produced by the con- 
tinual triggering of SF through compression in the gas 
peak (Lubow et al. 1986). Even if the spiral structure is 
not directly shock triggering SF (e.g., Foyle et al. 2010), 
provided there is an appreciable amount of SF in the 
arms, one still expects bright stars to be observed down- 
stream from the arms and this has been noted in several 
cases (e.g., Tamburro et al. 2008; Egusa et al. 2009). 

Roberts et al. (1975) found observational evidence for 
such offsets from a sample of 24 galaxies and posited that 
the shock strength and spiral strength could account for 
the ordering of spiral disks into Hubble type and lumi- 
nosity class. Figure [1] taken from R69 illustrates the po- 
sitions of the gas shock and young stars inside corotation 
of a trailing spiral pattern and illustrates the predicted 
offsets. Using a rigidly rotating spiral potential in hydro- 
dynamical simulations, Gittins & Clarke (2004) showed 
that one can find three spiral patterns each tracing a 
different stage in this process. The first is traced by 
the stellar mass density and the location of the poten- 
tial minima. The second spiral pattern is traced by the 
dust, which is known to be in locations of compressed 
gas and the third marks the location of the young stars. 
Due to the finite time for SF, the young stars should be 
downstream of the dense gas in this picture. The loca- 
tion of the compressed gas relative to the minimum is 
more complex; traditionally it is located upstream of the 
minimum, but the location of the shock depends on the 
sound speed of the gas and can be downstream of the 
minimum in a cold, or multiphase medium (Wada 2008; 
Dobbs 2007). 

Figure [2] shows what would be expected using a simple 
toy model of a two-arm spiral pattern with a corotation 
radius of 2.7Rexp (^-e. Kranz et al. 2003). We consider a 
gas and recent SF tracer with an onset time of 3 Myr be- 
tween the two. The upper right panel of Figure [2] shows 
the radial profile of the angular offsets between them 
calculated from Eq. [6l discussed in the following section, 
using the rotation curve shown on the left. In the inner 
regions large positive offsets are expected and beyond 
corotation the offsets become negative. The lower panel 
shows the position of the tracers in the spiral arms. The 
variation in the angular offsets of the tracers with radius, 
produces a spiral pattern that is effectively more tightly 
wound. The offsets are small enough that the arm pat- 
terns of the two tracers will likely overlap, depending on 
the amount of dispersion and width of the tracers. This 
could make by-eye determinations of offsets challenging. 



A number of observational studies have looked for qual- 
itative evidence for angular offsets of SF tracers. Early 
observations showed dust lanes on the inside part of the 
arms with respect to H II regions {e.g., Lynds 1970). 
Mathewson et al. (1972) and Rots (1975) found that 
atomic hydrogen was offset from optical arms and dust 
lanes in M51 and M81 respectively. A detailed study 
of M33 by Humphreys & Sandage (1980), however, pre- 
sented a more complex picture. While two of the arms 
did show dust lanes on the inner edge of the arm and 
bright young stars on the outer part of the arm, the third 
arm showed no such features or ordering. 

Ideally one would like to examine offsets between the 
gas and young stars. Due to the lack of atomic hydrogen 
(H I) in the inner part of the disk and the fact that the 
molecular gas forms giant molecular clouds out of which 
stars form, it is natural to look for offsets between CO 
and Ha or 24 /urn emission (with both Ha and 24 /am 
emission being tracers of recent SF, e.g., Calzetti et al. 
2007). A number of studies have looked for such offsets 
and examined the streaming motions near the arms, es- 
pecially in M51 and M81 {e.g., Vogel et al. 1988; Rand 
& Kulkarni 1990; Lord & Young 1990; Garcia-Burillo et 
al. 1993; Rand 1995; Loinard et al. 1996; Shetty et al. 
2007; Egusa et al. 2009). However, all such studies used 
by-eye determinations of these offsets by selecting indi- 
vidual patches of a single arm. Such methods can easily 
introduce potential biases. 

Tamburro et al. (2008, hereafter T08) developed an al- 
gorithmic technique to measure offsets between SF trac- 
ers by locating the peak of the cross-correlation function 
between H I and 24 jam emission. This also allowed for 
a measure of the timescale of going from H I atomic gas 
to dust-enshrouded massive stars (as traced by 24 /im 
emission). Their findings were in agreement with the 
simple prescription of R69. They found corotation radii 
at ~2.7rexp, which is consistent with other works {e.g., 
Kranz et al. 2003) and short timescales of 1-4 Myr for 
the timescales of SF emerging from H I. Indeed the an- 
gular offsets were found to be very small (5°) between 
the tracers. T08's work is the first to approach offset 
measurements using an algorithmic method as opposed 
to by-eye determinations and they used the highest res- 
olution images to date. Our method is based on the one 
developed by T08 and we describe it in greater detail in 
§2. 

Recently, Dobbs & Pringle (2010, hereafter DPIO) have 
looked in detail at the distribution of cluster ages across 
spiral arms in realistic simulations. In their sample, they 
included a galaxy with a fixed spiral potential which 
mimics a long-lived spiral with a constant pattern speed. 
The top left panel of Figure [3] shows the cluster age dis- 
tribution in this galaxy. If one compares this to the illus- 
tration of R69 (Figure [T]), one sees that the distribution is 
more complex in these simulations. However, the clusters 
are still ordered in the way expected. Previous studies 
by Dobbs et al. (2006) have also shown that the distri- 
bution of the giant molecular clouds has a sharp edge 
on the upstream side of the arms where fresh material 
is flowing in and a dispersed, smooth distribution on the 
downstream side. 

However, imposed stationary potentials are not realis- 
tic. DPIO also included three spiral structures without 
external potentials, in which the spiral structure forms 
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Figure 1. Illustration of the relative location of gas shock, sharp 
H I peak and newly formed stars in H II regions in the gaseous 
spiral arms of a trailing spiral pattern (top; Figure 7 of R69). 

naturally. They included a flocculent spiral as well as a 
barred and interacting galaxy. Unlike imposed stationary 
potentials, the active potentials do not produce a spatial 
separation between the gas and SF tracers. Dobbs & 
Bonnell (2008) have shown that the stars and gas are 
coincident with the spiral potential minima, because gas 
accumulates in the potential minima as the spiral struc- 
ture forms. The spiral arms in the flocculent model of 
DP 10 evolve over timescales of 100 Myr, which is typi- 
cal for spirals with recurrent spiral structure (Sellwood 
& Carlberg 1984; Fujii et al. 2011; Sellwood 2010). The 
spiral patterns found with active potentials have no fixed 
pattern speed or corotation radius. Longer lived tran- 
sient spirals, where spiral arms evolve over a few hun- 
dreds of Myr are possible {e.g. Thomasson et al. 1990), 
and may also lead to offsets similar to the static potential 
case. However such structures are not readily produced, 
at least in numerical simulations (Sellwood 2010), and we 
did not include this scenario in our simulated galaxies. 

In this work, we look for evidence for a systematic 
ordering of cold gas and young stars, which would be 
predicted by a long-lived spiral pattern with a constant 
pattern speed. We first examine the four simulations 
of DP 10 for angular offsets between the gas and stellar 
clusters and then examine a sample of observed galax- 
ies for offsets between various gas and SF tracers. In §2 
we describe the algorithmic method developed by T08 
to measure angular offsets between gas and SF tracers. 
In §3 we use this method to measure angular offsets in 
the simulations of DP 10, which cover a range of spiral 
structure formation mechanisms including a stationary 
spiral potential, bar potential, interacting galaxy and 
transient spiral structure. In §4 we turn to observations 
and present our sample of 12 nearby spiral galaxies with 
high-quality maps of cold gas and SF tracers including 
H I (THINGS), CO (HERACLES), 24 /im (SINGS), UV 
(GALEX) and 3.6 /im (IRAC). Using the same method 
employed for the simulations, we measure angular offsets 
between these observed tracers. In §5 we compare the 
observed cross-correlation functions and angular offsets 
to those in the simulations and assess the predictions of 
the model of a long-lived spiral structure with a constant 
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Figure 2. Toy model illustration showing the relative position 
between the gas and a recent SF tracer with an onset time of 3 
Myr between the two. The rotation curve of the galaxy is modeled 
with an arctan function (upper left) and the corotation of the spiral 
pattern is set at 2.7Rexp (dashed line). With a t^i? of 3 Myr, Eq.[6] 
determines the radial profile of the angular offsets between the two 
tracers (upper right). The relative positions of the spiral patterns 
of each tracer are shown below. 

pattern speed. In §6 we present our conclusions. 

2. MEASURING ANGULAR OFFSETS 

We seek to algorithmically measure angular offsets be- 
tween the gas and SF tracers in both simulations and 
observations, which are predicted by a long-lived, quasi- 
stationary, spiral structure theory. We use the method of 
Tamburro et al. 2008, which was developed to measure 
offsets between H I and 24 /im emission. Angular offsets 
between H I and 24 /im emission allow one to measure 
the timescale between gas compression and massive SF 
because peaks of H I column density are thought to be 
the sites of molecular cloud formation and the 24 /im 
emission traces the young dust enshrouded stars. We 
use this method to measure angular offsets between the 
gas and young stellar clusters in the simulations of DP 10 
and also extend T08's study to measure angular offsets 
between CO, which traces the molecular gas with 24 /im 
emission. This provides a more direct measure than that 
of H I and 24 /im emission. While T08 also used BIMA 
SONG CO maps (Heifer et al. 2003) to measure angular 
offsets between the molecular gas and 24 /im, the low sen- 
sitivity of these images provided too few points to make 
accurate estimates of the timescale. We also include the 
cross-correlation of the gas tracers with UV images and 
3.6 /im images which probe longer timescales. We first 
describe the method of T08 in detail. 

2.1. Method of Tamburro et al. 2008 

T08 measured angular offsets by locating the peak of 
the cross-correlation function between H I and 24 /im 
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emission in radial annuli for 14 disk galaxies. 

The time between these two events, tHi^24fim^ the lo- 
cal circular velocity, Vc{r), and a pattern speed, l^p, for 
the spiral arms, lead to an angular offset between these 
two tracers of: 



A(/)(r) = {Q{r) - %)tHi- 



► 24/im- 



(1) 



Thus, this method is independent of the number of 
arms and shape of the spiral pattern, but it does assume 
a fixed pattern speed. The angular offsets should vary 
with radius. When the disk rotates faster than the pat- 
tern (inside corotation) we expect A(/) > and when the 
pattern rotates faster than the disk (beyond corotation) 
we expect < 0. At corotation Vt{Rcor) = ^p, the sign 
changes and A(j) = 0. 

In order to measure A(/), T08 divided the images into 
radial annuli and cross-correlated them to find the rel- 
ative offset with the best match. The best match be- 
tween the two images is found by minimizing the follow- 
ing quantity as a function of the shift, / (lag): 



(0 = ^k[xk - Vk-i] 



(2) 



The sum is carried out over all N elements of x, where 
X and y refer to the two tracers at a given radius as a 
function of azimuth, (/), such that: 

Xk = fHi{(t>k\r) and yu-i = f2Aii7n{(t>k-i\r). (3) 

By minimizing one maximizes the following quantity: 



CCx,y{l) = ^k[Xk X Vk-l] 



(4) 



defined as the cross-correlation coefficient. T08 used the 
normalized version such that: 



T.k[{xk - x){yk- 



y^T^kixk - xYY^kiyk - vY 



(5) 



where x and y are the means of x and y respectively. 
In this way, perfectly identical patterns have a cross- 
correlation coefficient of unity and highly dissimilar pat- 
terns have a value much less than one. At each an- 
nulus, T08 searched for a local maximum of the cc(l) 
around / :^ 0. The location of the maximum, Imaxi de- 
fines the angular offset between the two tracers such that 
A0(r) = —Imaxij)' If the value of the cc{l) peak was less 
than 0.3, the point was rejected. 

With the angular offset, A0(r), and the rotation curve, 
Vc{r)^ one can write Eq. [T]as: 



A^(r) 



Vc{r) 



\tHI- 



►24/im- 



(6) 



T08 used fitting of the above with the measured an- 
gular offsets and derived best fits for tHi^2A^im and Vtp. 
This method can be used for any two tracers. We first 
employ this method on the simulations of DP 10 and ex- 
amine which models fit the predictions of R69. We then 
use this method to repeat the work of T08 as well as look 
for offsets between 24 jiia and the higher sensitivity CO 
imaging from HERACLES, which traces the molecular 
gas. We also examine UV and 3.6 jiia for offsets with 
respect to the molecular gas. 




Figure 3. Distribution of stellar cluster ages (shown in different 
colors) in the spiral arms of a long-lived spiral potential with a 
constant pattern speed (top left), a barred galaxy (top right), a self- 
excited flocculent galaxy (bottom left), and a galaxy interacting 
with a companion (bottom right; Figure 2 of DP 10). While the 
ordering of the stellar clusters is systematic in the long-lived spiral, 
the other cases show a much more complex distribution. 



3. ANGULAR OFFSETS IN SIMULATIONS 

DP 10 have recently simulated four different models of 
spiral structure and compared the distribution of stel- 
lar clusters of different ages. Their simulations included 
a stationary spiral potential which mimics a long-lived 
spiral structure with a constant pattern speed, a barred 
galaxy, an interacting galaxy like NGC 5194 and a tran- 
sient spiral produced using the spiral potential from the 
models of Sellwood & Carlberg (1984). As we saw in 
§1, the long-lived spiral structure was simulated using 
a stationary potential and indeed the ordering of the 
stellar cluster ages fit the predictions of R69 and others 
(see Figure [1]). Since the long-lived structure requires 
a stationary potential in order to ensure its longevity, 
the other three cases may be more realistic (see Figure [3] 
for all four cases). In these cases the ordering is much 
more complex. For the fiocculent galaxy, each different 
segment of the spiral arms contains clusters of different 
ages. The interacting galaxy has an almost incoherent 
distribution. 

We have used the same simulations presented in DP 10 
and analyzed them for angular offsets between the gas 
and stellar clusters of different ages in the way designed 
by T08 for observations. For different snapshots we con- 
vert the gas and stellar cluster onto an R-^ grid and cross 
correlate the annuli and locate the peak of the cross- 
correlation function as described in §2. We show here 
the results of the cross-correlation between the gas and 
100 Myr clusters in all four simulations. The gas and 
the 100 Myr clusters represent timescales close to those 
observed between H2 and the UV, which we will examine 
in §5. 

3.1. Stationary Spiral Potential 
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Figure 4. Cross correlation function between the gas and the 
100 Myr clusters for the simulated stationary spiral potential (top 
panel), bar potential (second from top), interacting spiral (third 
from top) and flocculent spiral (bottom) at a flducial radius of 7 
kpc in each case. The cross-correlation function of the stationary 
potential shows four equally spaced peaks due to the self-similarity 
of the four-armed spiral structure. In the other cases, the central 
peaks are less well-deflned and show very small, if any, offsets. 

The stationary spiral potential is an imposed potential 
with a constant pattern speed. Thus, we expect an or- 
dering of the clusters and gas to be in the way predicted. 
The pattern has a corotation radius at 11 kpc. The up- 
per panel Figure H] shows the cross-correlation function 
between the gas and the expected distribution of the 100 
Myr stellar clusters at a fiducial radius of 7 kpc. The 
galaxy has four arms, so there are four peaks in the cross- 
correlation function due to the self-similarity of the pat- 
tern. At each radius, we fit the central peak with a four 
degree polynomial, as was done by T08. The peak of the 
polynomial locates the angular offset. Figure [5] shows 
the chosen angular offset between the gas and 100 Myr 
clusters at a series of radial annuli. As one moves out in 
radius the offsets decrease in value as we saw with the toy 
model in §1. Unfortunately, the stellar clusters do not ex- 
tend beyond corotation so we are unable to confirm that 
the sign of the offsets changes beyond corotation. 

Since the arms are well-defined in the stationary po- 
tential we fit straight lines to the gas arms in logR-(^ and 
defined a ridge line to the arm. At each radial position 
the arm is at a position (j)ridge' We then measure the dis- 
tance of the clusters with respect to the arm. Figure [6] 




Figure 5. Radial proflle of the angular offsets based on the peak of 
the cross-correlation function between the gas and 100 Myr cluster 
particles for the stationary spiral potential. The offsets show a 
systematic trend going from high to low values. The gas and star 
particles in this simulation do not extend beyond the corotation of 
the spiral potential. Thus, we are unable to currently test whether 
the offsets cross zero at this point. 

shows stellar clusters at 2, 50, and 100 Myr and their an- 
gular distance to the ridge line. We have only fit one arm 
so the self-repeating pattern is due to the other arms. 
One notes that the older clusters have drifted further 
downstream from the arms {i.e. {(j)ridge — ^) increases 
for clusters of older ages). Both these findings and the 
systematic ordering of the angular offsets fit the predic- 
tions of a long-lived, stable, spiral structure. This is not 
surprising because in these simulations such a structure 
is imposed. This exercise shows that if such a structure 
exists in nature, we should be able to detect offsets be- 
tween gas and SF tracers using our analysis. 

We now turn to the other simulations, where the spiral 
structure has formed naturally in the simulations and is 
not fixed with a stationary potential. In these cases, 
however, we do not expect a constant pattern speed for 
the spiral pattern, which is one of the assumptions in the 
R69 model. 

3.2. Barred Galaxy 

The second panel from the top of Figure |4] shows the 
cross-correlation functions of the gas and the expected 
distribution of the 100 Myr stellar clusters for the barred 
potential at a fiducial radius of 7 kpc. Here we find 
two cross-correlation peaks due to the two arms in this 
galaxy, but the cross-correlation peaks are small and, in 
many cases, are below the 0.3 cutoff value used by T08 
and which we also adopt for the observed galaxies. De- 
spite these low cc{l) values, the peak is well-defined so 
we are able to determine the offset angles. We find that 
the offsets are very small and, at all radii, are less than a 
degree. One concern is that by 100 Myr the ordering of 
the particles has been disrupted. In comparison to the 
other simulated galaxies, the bar primarily extends over 
smaller radii, and the orbital times are shorter. By 100 
Myr, the particles have made at least one and possibly 
two passages around the galaxy, disrupting the order- 
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Figure 6. Angular distance of clusters of varying ages to the 
gaseous spiral arms in the stationary spiral potential simulation. 
By fitting a ridge line to a spiral arm in the gas ((f)ridge)^ we mea- 
sure the angular distance to each age of stellar cluster (2 Myr: 
black; 50 Myr: blue; 100 Myr: red). As the clusters age they drift 
further downstream with respect to the arm, as expected in the 
stationary, spiral potential model. 

ing of the gas and star particles. Due to this, we also 
examined star particles at early times, before complete 
passages. For the 50 Myr clusters the offsets were easier 
to detect and larger, but still quite small with the largest 
being ~ 5°. A further challenge in the barred galaxy is 
the highly elliptical orbits of the particles, which make 
the transitions more difficult to detect using the cross- 
correlation method, which relies on the azimuthal angle 
of the particles. In general, we find that offsets can be 
detected for the barred galaxy, but that they are quite 
small and the ordering of the particles can be easily dis- 
rupted. 

3.3. Interacting Galaxy 

The second panel from the bottom of Figure H] shows 
the cross-correlation functions for the gas and the ex- 
pected distribution of the 100 Myr stellar clusters of the 
interacting galaxy, which is very similar to NGC 5194, 
at a fiducial radius of 7 kpc. In this case, we find only a 
single, broad, uneven cross-correlation peak. The mea- 
sured angular offsets show no trend and in most cases 
the angular offsets are close to zero and thus we do not 
show them graphically. Since this simulated galaxy is 
very similar to NGC 5194 (see DP 10 for further details), 
it is interesting to compare the cross-correlation function 
with the observations, which we will do in the following 
section. 

3.4. Transient, Flocculent Spiral 

Finally we also examine the cross-correlation function 
of the transient, more fiocculent spiral structure that 



used the spiral potential from the models of Sellwood 
& Carlberg (1984). The bottom panel of Figure |4] shows 
the cross-correlation function between the gas and the 
expected distribution of the 100 Myr stellar clusters at a 
fiducial radius of 7 kpc. At this radius and all other ones 
no obvious cross-correlation peak is found, suggesting no 
preferred position between the gas and clusters. We were 
unable to properly fit a peak and measure angular offsets 
and thus, we do not show the radial variation here. 

3.5. Summary of Numerical Models 

If a long-lived spiral structure with a constant pat- 
tern speed accurately describes spiral structure, than we 
should be able to detect angular offsets between gas and 
SF tracers as we have seen with the imposed stationary 
potential of DP 10. However, if a stationary potential is 
not imposed, simulations of spiral structure do not show 
a systematic ordering of angular offsets between the gas 
and stellar clusters in fitting with the model predictions 
of a long-lived, stable, spiral structure. In the case of 
the barred galaxy, the cross-correlation peak could be 
fit, but the values were very low and the angular offsets 
measured, were small. In the other two cases, the offsets 
showed no radial variation and were close to zero. We 
now turn to our sample of observed galaxies and compare 
our findings. 

4. ANGULAR OFFSETS IN OBSERVED GALAXIES 

We chose 12 galaxies in common with the sample of 
T08 that have coverage in THINGS (Walter et al. 2008) 
and SINGS (Kennicutt et al. 2003). Eight of these galax- 
ies also have coverage in HERACLES (Leroy et al. 2009). 
For three galaxies we also included FUV images from 
GALEX (Gil de Paz et al. 2007) and 3.6 /im images from 
IRAC (Kennicutt et al. 2003). The sample presents a 
mixture of both grand design and more fiocculent spi- 
rals, which are similar to the four simulated galaxies of 
DPIO. 

The galaxies were deprojected according to the values 
of T08, with the exception of NGC 2841, NGC 3521 and 
NGC 5194 (see Table 1). For NGC 5194, we used the val- 
ues adopted by Leroy et al. (2008). NGC 2841 and NGC 
3521 were deprojected using GALFIT (Peng et al. 2002) 
and are roughly equivalent to those found in HyperLeda 
Extragalactic Database. Figure [71 shows deprojected 3.6 
/im images for our sample as well as the direction of 
rotation, which is defined by assuming a trailing spiral 
pattern. 

The images were aligned to the THINGS astrometric 
grid and were degraded to a common resolution of either 
6'' for our analysis of H I and 24 /im images or 13'' for our 
analysis of CO, 24 /im, UV and 3.6 /im images. Foyle et 
al. (2010) describes in detail the processing of the images. 

After initial processing, the images are translated from 
rectangular to polar coordinates, (r, ^), and divided into 
annuli of 5" width as was done by T08. At our convolved 
resolution, this allows us to Nyquist sample the data. As 
an example. Figure [8] shows the H I emission (bottom 
panel) and H2 map (upper panel) in polar coordinates for 
NGC 5194 with the contours of 24 /im emission overlaid. 
One notes that in both cases the arm patterns in these 
tracers are nearly coincident. Thus, any offsets between 
these tracers must be small. 
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Table 1 

Properties of 12 Sample Galaxies 



Name 


Inclination 


P.A. 


D 


ri 






n 


n 


[Mpc] 


[km s"-*^ 


NGC 628 


7 


20 


7.3 


1.1 


220 


NGC 2403 


63 


124 


3.22 


1.3 


128 


NGC 2841 


63 


148 


14.1 


0.92 


331 


NGC 3031 


59 


330 


3.63 


3.62 


256 


NGC 3351 


41 


192 


9.33 


0.86 


210 


NGC 3521 


65 


162 


10.05 


0.74 


242 


NGC 3621 


65 


345 


6.64 


0.8 


144 


NGC 3627 


62 


173 


9.25 


0.95 


204 


NGC 5055 


59 


102 


7.82 


1.16 


209 


NGC 5194 


20 


172 


7.77 


1.39 


242 


NGC 6946 


33 


242 


5.5 


1.73 


201 


NGC 7793 


50 


290 


3.82 


1.16 


109 



Note. — Sample properties: the inclination and position an- 
gles used to deproject the galaxies, their adopted distance, scale 
length and maximum amplitude of their rotation velocity. With 
the exception of NGC 2841, NGC 3521 and NGC 5194, all values 
are those from T08. 
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5055 
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3351 
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3621 



2403 



2841 



3031 




Figure 7. Sample of 12 galaxies (numbers correspond to NGC 
numbers). We show the deprojected 3.6 /xm images with arrows 
denoting the direction of rotation based on the assumption of a 
trailing spiral pattern. 



The cross-correlation coefficients were calculated for 
angular offsets from -180° to +180° in increments of ~ 
0.1°. We restrict the search for the peak between ±30°, 
because the offsets must be at least this small based on 
visual inspection (see Figure [8j). We fit the region of 
the maxima with a four degree polynomial and calculate 
its peak value. The angular position of the peak is the 
angular offset, A^(r), between the tracers at each radius. 
Any local maxima with a cross-correlation coefficient less 
than 0.3 is rejected, in following T08. The sign of the 
angular offsets is corrected for the direction of rotation 
of the galaxy. 

4.1. Angular Offsets Between HI and 24 l-im Images 

Figure [9] shows a sample cross-correlation function be- 
tween H I and 24 /im at a radial annulus of 80^' for NGC 
5194. One notes two broad peaks due to the two spiral 
arms and the self-similarity of the structure. The zoomed 
box on the right shows the peak and the polynomial fit 
(gray curve), the maximum of which locates the angular 
offset at that radius (vertical black line). The offsets are 
measured in a similar fashion at all radii. 

Figure [10] shows the angular offsets (black points) ver- 
sus radius for all 12 galaxies. The colors reflect the value 
of the cross-correlation with low values (^ 0.3) being 
blue and high values 0.8) being red in the neighbor- 
hood of the chosen offset corresponding to the peak value 
(±20°). Around the peak value, there is a broad region 
where cc{l) is high. One expects to find positive angular 
offsets smoothly decreasing with radius and crossing zero 
at corotation. The radial profiles of the angular offsets do 
not show any smooth trend of going from positive to neg- 
ative values and show considerable scatter. The profiles 
do not correspond in any way to the model predictions. 
Our profiles also do not agree with the profiles found by 
T08, who found offsets that agreed with the picture of 
R69. In §4.3 we examine this in detail. However, if the 
reader is not interested in this discussion, the §4.3 can 
be skipped. 

4.2. Angular Offsets Between II2 and 24 IJim, UV and 

3.6 jam Images 

The H I distribution is not as concentrated to the spiral 
arms as H2 {e.g., Foyle et al. 2010). Since the H I emis- 
sion is more evenly distributed than the 24 jam emission, 
this may make the measurements of angular offsets more 
difficult. Thus, we cross-correlate H2 maps with other 
SF tracers to measure angular offsets. 

Eight of the galaxies had CO coverage from HERA- 
CLES at 13'' resolution. As described in Foyle et al. 
(2010), we process and convert those maps to produce 
deprojected maps of H2 using an Xco conversion factor. 
We then transform the images to polar coordinates. As 
we did with H I and 24 /im emission, we measure angu- 
lar offsets between the H2 and 24 /im emission by locat- 
ing the peak of the cross-correlation function. Figure [TTI 
shows the radial profile of the angular offsets between 
these two tracers. In a similar fashion to Figure [TOl we 
find little evidence for systematic offsets, in contrast to 
R69 and T08. There is some evidence for a sequence of 
offsets in NGC 6946 and NGC 5194. Thus, we examine 
these galaxies as well as NGC 628 in greater detail and 
with other tracers for possible angular offsets. 



8 



Foyle et al. 




Figure 8. Hi (bottom) and CO (top) emission of NGC 5194 in polar coordinates with contours of the 24 /im emission overlaid. It is 
clear from these overlays that any angular offsets between both sets of tracers must be small. 



NGC 5194 

■ ■ ] 0.8 




-180 -90 90 180 -30 30 
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Figure 9. Example of a cross-correlation function, cc{l), for H I 
and 24 /xm as a function of azimuthal offset for NGC 5194 at a ra- 
dius of 80" . Two broad peaks are present due to the self-similarity 
of the two spiral arms. The box on the right shows the region 
zoomed in (±30°). The fourth degree polynomial fit overlaid in 
gray. The location of the maximum of the polynomial (black ver- 
tical line) gives the angular offset between the two tracers. There 
is negligible offset between these tracers at this radius (2.7°). This 
example is not atypical for different galaxies, and radii. 

For three galaxies with prominent spiral structure we 
added FUV images from GALEX (Gil de Paz et al. 2007) 
which trace a longer timescale and should show the great- 
est offsets. We also included 3.6 jiia emission maps which 
trace the underlying old stellar population. We mea- 
sured angular offsets between the H2 and the UV emis- 
sion as well as H2 and the 3.6 jim. emission using the 
cross-correlation method described above. The timescale 
between the H2 and the UV is similar to the timescale 
between the gas and the 100 Myr clusters in the simula- 
tions of DPIO (see §3). 



Figure [12] shows the angular offset profiles for cross- 
correlations of H2 with the UV emission on the left and 
the H2 with the 3.6 /im emission on the right. Of these 
galaxies, only NGC 628 presents any trend suggestive 
of the picture of R69. However, it is important to note 
that many of these offsets are well below the resolution 
of the images, especially for the case of H2, which only 
has 13'' resolution. We also tried other combinations of 
these images and divided the images into quadrants to 
isolate individual arms. Our findings were similar in all 
cases. 

4.3. Comparison with Tamburro et al. (2008) 

We have found no evidence for a systematic ordering 
of angular offsets in the way found by T08. We examine 
the cause of this discrepancy in greater detail here. Fig- 
ure [13] shows the cross-correlation functions of H I with 
24 /im emission at several radial annuli for NGC 5194. 
The left panels show the fuh range from -180° to +180°. 
We also zoom in on the peak of the function in the boxed 
regions and display these on the right. The zoomed re- 
gions reflect the range of -30° to +30°. The grey curve 
shows the fourth degree polynomial fit. 

The top panel of Figure [13] exhibits what one would 
hope to find when seeking the maximum of cc{l). At this 
radius a single, narrow peak is found. The other panels, 
show, however, that at most radii, the cc{l) is more com- 
plex, with a very broad peak and uneven features. The 
polynomial fit is sensitive to the range over which one fits 
and selecting this range can be subjective. Due to this, 
we also tried selecting the maximum cc{l) value, rather 
than fitting the peak. Even in this case, no systematic 
variation of offsets were found. We tried a number of 
variations of T08's method, including smoothing and fil- 
tering of the images, variations of the cross-correlation 
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Figure 10. Radial profiles of the angular offsets based on the 
location of the peak of the cross-correlation function between H I 
and 24 /im vs. radius (black dots). The colors reflect the value 
of the cross-correlation coefficients from high (red) to low (blue) 
in the range of ±20°. Values below 0.3, our chosen cutoff, are 
not shown. The angular offsets are corrected for the direction of 
rotation. None of the galaxies show the trend predicted by the R69 
model. 

calculation and restrictions of narrow regions around the 
arms in the case of non-axisymmetric distributions of 
stars and gas and were unable to reproduce the findings 
of T08. 

We carried out a careful comparison of NGC 628 also 
studied by T08. While we were able to reproduce some 
of their offsets, the offsets measured were sensitive to 
small changes in the polynomial fit and the way in which 
the cross-correlation function is calculated {i.e. use of 
mean or median for x and y in Eq. [5j) . Figure [M] shows 
the offsets found by T08 (black) and two of our attempts 
using different fitting ranges for the polynomial (blue and 
red). While T08's offsets show a clear progression from 
high values to low values, which cross zero, our points 
are much more scattered. In part, we can attribute this 
to the fact that T08 selected the peak to fit and chose 
the range over which they fit by-eye. We also tried to 
choose specific fitting ranges at each radius, but found 
the choice to be subjective. 

Beyond the sensitivity of the peak position to the range 
of the polynomial fit, resolution may be a cause for the 
lack of agreement with T08. The angular resolution 
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Figure 11. Same as Figure [Tol but showing the location of the 
peak of the cross-correlation function between H2 and 24 /im emis- 
sion. 

varies with radius and the thin black line in Figure [HI 
shows that all of the offsets found by T08 are below the 
resolution of the image. The common image resolution 
is 6'' and at a radius of 30^' this corresponds to an an- 
gle of 11°. Most of the offsets found by T08 were less 
than 5° except in the very inner regions, where resolu- 
tion concerns are even greater. However, given that the 
cross-correlation function is calculated considering all az- 
imuthal positions, it is possible that one can probe below 
the resolution limit to measure the offsets. The fact that 
T08's offsets are not randomly scattered points to this 
possibility. 

It is not clear how resolution might affect the cross- 
correlation function and the determination of the peak. 
We first examine how resolution affects the autocorrela- 
tion function of two H I maps at a resolution of 6'' and 
13'' using NGC 5194. The autocorrelation function is 
calculated like the cc{l) (Eq. [5j), but in this case both x 
and y represent the same H I map. For both resolutions, 
we expect the peak of the autocorrelation function to be 
centered at zero with a maximum of unity. Figure [T5l 
shows the autocorrelation function for a series of annuli. 
Since the azimuthal resolution varies with radius, we list 
the azimuthal resolution for each pair of images in the 
upper left of the panels. We see that, as expected, the 
autocorrelation functions have a peak at zero and max- 
ima of unity. However, one notes that as the resolution 
size increases, the autocorrelation function broadens and 
has a higher coefficient over a greater range of angles. We 
divided each annulus into units smaller than the resolu- 
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Figure 12. Same as Figure fTol but showing the location of the 
peak of the cross-correlation function between H2 and UV emission 
(left) and H2 and 3.6 /xm emission (right). 
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Figure 13. Cross-correlation functions cc(l) between H I and 24 
/xm as a function of azimuth offset for NGC 5194. Each panel 
represents neighboring annuli from 60^^ (bottom) to 80'^ (top). The 
right column shows the boxed region zoomed in (±30°). The fourth 
degree polynomial fits are overlaid in gray. The peak is often quite 
broad and bumpy and the polynomial peak depends heavily on the 
range of points selected. 



tion. Thus, as the resolution size increases, the fraction 
of pixels with identical values increases. This provides a 
better agreement in the cross-correlation, which broad- 
ens and increases the peak. 

In order to examine how the peak position might vary 
with resolution, we consider the cross-correlation of H I 
and 24 /am at two different resolutions in Figure [161 
We use pairs of images at our common resolution of 6" 
(black), and 13^' (gray). Like the autocorrelation func- 
tion, as the resolution size increases, the cc(/) broadens 
and has a higher value. The functions are very similar 
in shape, but as one can see in the zoomed in regions, 
when the resolution size is large (small radii), the shape 
of the peak can change substantially. This affects the 
polynomial fit and ultimately the position of the peak 
center. In the inner regions, where the resolution effects 
are greatest, we found that the peak positions varied by 
as much as 8°. However, at larger radii the difference was 
less than a degree and the median deviation was 0.5°. 
Thus, while the peak position remains roughly the same, 
it is important to consider the uncertainties introduced 
by resolution in determining its precise position. 

In summary, we can attribute our lack of agreement 
with T08 to a combination of three causes: 1) the offset 
value is very sensitive to the polynomial fitting of the 
peak; 2) the offset value is very sensitive to small changes 
in the way the cc{l) function is calculated; and 3) almost 
all of the offsets are below the resolution limit of the 
images and we have found that changes in the resolution 
of the image do affect the shape of the cc{l) function 
particularly at small radii. 
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Figure 14. Comparison of the offsets found by T08 (black) and 
the offsets measured in this study using two slightly different fitting 
ranges of the cc(l) (blue and red). The offsets measured by T08 
are all below the angular resolution limit, which varies with radius 
(thin black line). 

5. COMPARISON OF SIMULATIONS AND OBSERVATIONS 

For all galaxies in the observed sample and between 
all tracers we found a clear peak in the cross-correlation 
function (see Fig [9]). The location of the peak marks 
the angular offset between the two tracers considered. 
However, the measured angular offsets showed no sys- 
tematic variation as would be expected from a long-lived, 
quasi-stationary spiral structure. Among the simulated 
galaxies only the one with an imposed stationary po- 
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Figure 15. Effect of spatial resolution on the autocorrelation 
function of H I maps with 6^' (black) and 13^^ (gray) resolution 
for NGC 5194 at a series of fiducial radii. Since the azimuthal res- 
olution varies with radius, it is listed in the upper left of each panel. 
As one moves out in radius (from bo ttom to top) the azimuthal 
resolution size increases. As in Figure [T3] we zoom in on the peak 
between zb 30°. We see that as the size of azimuthal resolution 
increases, the peak of cc(l) broadens and increases in value. 



tential showed clear offsets with a systematic variation. 
The barred, interacting and flocculent transient spirals 
all showed very small offsets (less than a few degrees) and 
the interacting and flocculent spirals showed no clear sys- 
tematic variation. This agrees well with what we found in 
the observations for the various tracers of gas and recent 
SF. The absence of observed offsets not only suggests 
that these galaxies do not have a static spiral pattern, 
but also that their structure is unlikely to be the result 
of a single mode transient spiral that persists for longer 
times. 

However, upon comparing the shape of the cross- 
correlation functions of the simulated galaxies, with 
those observed, we do not find a simulated model that 
agrees in detail. The observations tend to show at least 
two broad peaks, depending on the nature and number 
of arms in the spiral pattern (see Fig [9j) . The simulated 
barred galaxy, did show two clear peaks, but the cross- 
correlation coefficients were much smaller than those seen 
in observations. In the case of the flocculent and interact- 
ing spiral simulations, we did not find such clear features. 
It is interesting to compare the cross-correlation function 
shown in Fig[9]with that of the interacting, simulated spi- 
ral galaxy in Fig IH The interacting spiral is structurally 
similar to NGC 5194. While the cross-correlation func- 
tion of NGC 5194 showed two broad peaks, the simulated 
galaxy showed only one. One possible cause of the dis- 
crepancy in the shape cross-correlation function in the 
simulations and observations is because we only have a 
relatively small number of discrete clusters in the simula- 
tions, rather than the more continuous map of emission 
in the observation. 
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Figure 16. Effect of spatial resolution on the cc{l) function cal- 
culated using H I and 24 /im maps with 6^^ (black) and 13^^ (gray) 
resolution for NGC 5194 at a series of fiducial radii. Since the az- 
imuthal resolution varies with radius, it is listed in the upper left 
of each panel. As one moves out in radius (from bo ttom to top) 
the azimuthal resolution size increases. As in Figure fT3l we zoom 
in on the peak between ± 30°. While the cc{l) is roughly similar, 
features on the peak change substantially, which may introduce 
differences in the polynomial fit and the position of the maxima. 



6. CONCLUSIONS 

We have used simulations from DP 10 and high res- 
olution observations for a sample of 12 spiral galaxies 
to look for evidence for a model of spiral structure that 
is long-lived, quasi-stationary with a constant pattern 
speed. This model predicts that SF tracers should be 
offset from one another reflecting the relative velocity 
between the disk and pattern and the onset time for SF. 

We have used an algorithmic technique developed by 
T08 to measure angular offsets between the gas and stel- 
lar clusters in four simulated galaxies with different spi- 
ral structures including a galaxy with a stationary spiral 
potential, a barred galaxy, an interacting galaxy and a 
more flocculent, transient spiral structure. Only when a 
stationary spiral potential was imposed did we find clear 
angular offsets ordered in the way predicted by a long- 
lived, stable structure. For the barred galaxy, though we 
found angular offsets for younger clusters, and some in- 
dication of an age transition, the offsets were very small 
(<5°), and thus inconclusive. For the interacting and 
transient spirals, there were no clear cross correlation 
peaks and no measurable offsets. 

We then used the same method to measure angular off- 
sets between the gas and SF tracers in a sample of 12 ob- 
served galaxies. In all cases we did not have a systematic 
ordering in the way predicted by the simple prescription 
of R69. We specifically measured offsets between H I and 
24 jam emission in order to directly compare with the 
work of T08 and are unable to reproduce their results 
(see §4.3 for a discussion). This study has also measured 
offsets between CO, which traces the molecular gas with 
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the 24 /im emission. Even in this case, which directly 
probes the time between molecular cloud formation and 
SF, we find no systematic ordering of angular offsets. For 
three galaxies we also examined the cross-correlation of 
other tracer combinations including UV and 3.6 /im and 
found similar results. 

Our results contradict previous by-eye estimates of 
such offsets. While there may be some patches of SF 
tracers that show by-eye offsets, an algorithmic analy- 
sis shows that such patches are isolated and there is no 
overall trend which supports the prediction of a model 
where the spiral structure is long-lived with a constant 
pattern speed. Our technique made use of the highest 
quality images to date and, provided the timescales of 
SF are at least a few Myr, angular offsets should have 
been detectable. There are number of possible reasons 
why such a systematic ordering was not detected. It 
may be the case that there are no angular offsets be- 
tween the tracers. This would be the case if the spiral 
structure were simply a result of sheared patches of star 
forming regions or if the structure was rapidly dissolving 
and reforming. However, it is also possible that offsets 
are present, but that structure is more complex than a 
single pattern with a constant pattern speed. Indeed, re- 
cent studies have measured multiple pattern speeds for 
several spiral galaxies (Meidt et al. 2009). Furthermore, 
there is substantial emission from SF tracers in the in- 
terarm regions (at least 30%) (Foyle et al. 2010). This 
may be obfuscating the offsets near the spiral arms due 
to our cross- correlation method which uses all azimuthal 
positions. The arms may also have slightly different off- 
sets and the stars may have highly elliptical orbits, both 
of which would introduce uncertainties when fitting over 
the whole galaxy as our method relies on a polar coordi- 
nate system. Thus, there is still room for the possibility 
that systematic offsets between SF tracers exists in the 
arms, but we may not be able to quantitatively measure 
them with this technique. Furthermore, the fact a cross- 
correlation peak could be found at most radii, suggests 
that angular offset measurements could still be used to 
understand star-forming sequences if the local dynamics 
are known. 

The fact that there is little evidence for angular offsets 
between SF tracers as predicted by the model of a long- 
lived spiral structure with a constant pattern speed lends 
support to a model with a transient spiral structure that 
reorganizes the interstellar medium {e.g., Elmegreen & 
Elmegreen 1986; Dobbs & Bonneh 2008; Sellwood 2010) 
or models that have spiral arms forming due to the shear- 
ing of gas and stars by differential rotation {e.g., Sei- 
den & Gerola 1979 and Elmegreen et al. 2003). Sim- 
ulations largely support such pictures and other obser- 
vational studies have shown that most spiral structures 
are quite complex with multiple pattern speeds {e.g., 
Meidt et al. 2009). We caution, however, that the cross- 
correlation functions of the other simulated spiral struc- 
tures in this study, did not agree in detail with the obser- 
vations. Thus, continued detailed comparisons between 
observations and simulations will be required in order to 
uncover the nature and persistence of spiral structure. 
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